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Abstract 

We consider the orbital angular momentum of a free electron vortex moving 
in a uniform magnetic field. We identify three contributions to this angular mo¬ 
mentum: the canonical orbital angular momentum associated with the vortex, 
the angular momentum of the cyclotron orbit of the wavefunction, and a dia¬ 
magnetic angular momentum. The cyclotron and diamagnetic angular momenta 
are found to be separable according to the parallel axis theorem. This means 
that rotations can occur with respect to two or more axes simultaneously, which 
can be observed with superpositions of vortex states. 
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1 Introduction 

In classical mechanics, an electron moving in a uniform magnetic field follows a circular 
orbit in the plane perpendicular to the field, as dictated by the Lorentz force. This so- 
called cyclotron motion, which occurs even in the absence of a central potential, plays 
a role in areas as diverse as particle physics [1, 2], electron microscopy [3], plasma 
physics [4], and also in microwave ovens [5]. 

Electrons can also possess quantized canonical orbital angular momentum which 
does not depend on the presence of a magnetic field. This is well known in the case of 
bound states in atoms and quantum dots [6], in which there is a confining potential, 
however recently it was discovered that free electrons can be imprinted with orbital 
angular momentum. Electron vortex beams [7, 8, 9, 10], generated, for example, 
in electron microscopes, have twisted wavefronts, and resemble freely propagating 
atomic orbitals. The understanding, manipulation and exploitation of this angular 
momentum for a range of technological applications is currently a very active area of 
investigation [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. 

In a magnetic field, an electron can possess both canonical orbital angular mo¬ 
mentum, and angular momentum arising from the interaction with the field. If the 
magnetic field is uniform, the canonical angular momentum in the direction of the 
field is independent of the field and is constant [22]. Meanwhile, the magnetic field 
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induces an additional current within the electron’s wavefunction which gives rise to 
a diamagnetic angular momentum [23, 24, 25, 26]. Manipulating the canonical and 
diamagnetic orbital angular momenta of free electrons recently led to the first di¬ 
rect imaging of Landau states [27, 28]. In addition to the diamagnetic rotation of 
an electron’s wavefunction, however, in a magnetic field there will generally also be 
a cyclotron orbit of the centre of mass of the wavefunction [22, 29]. The angular 
momentum associated with this cyclotron motion has not previously been considered. 

In this paper, we show that the total orbital angular momentum of the electron 
is described by the parallel axis theorem. This angular momentum comprises the 
canonical and diamagnetic components, which are associated with rotation relative 
to the centre of mass of the wavefunction, and a cyclotron component which has 
expectation value equal to that for the classical orbit. Interestingly, for free electrons 
all three of these components can have similar magnitude. This means that the 
trajectory of the electron is strongly dependent on how these angular momenta add 
and subtract. Further, we show that different cyclotron orbits can be superposed, 
leading to rotations with respect to multiple parallel axes, and periodic interference. 
Our results suggest novel means of structuring electron beams for use in specific 
applications, such as probing magnetic and chiral materials. 


2 Model 


We consider an otherwise free electron moving under the influence of a uniform mag¬ 
netic field. We take the direction of this magnetic field to define our 2 : axis, and 
consider the motion of the electron within the x-y plane. The electron may also be 
moving in the z direction, however the component of its momentum in this direction is 
a constant of motion [25], and will not affect our results. We consider non-relativistic 
energies, meaning that the spin angular momentum is also constant, and can be sep¬ 
arated from the orbital motion of the electron [30]. In what follows, we shall consider 
only the electron’s orbital angular momentum. The magnetic field B = Bz can be 
described by the cylindrically symmetric vector potential A = Bp^/2, where p and 0 
are cylindrical polar coordinates. This choice of gauge is convenient as it allows us to 
exploit the rotational symmetry of the magnetic field. Our Hamiltonian is therefore 


H = 


{pT? 

2m 


1 

2m 


{pT? + - eBLl 


( 1 ) 


where = mv± = — is the component of the electron’s kinetic momentum 

in the plane perpendicular to the magnetic field, p^f^ = —ih'V± is the corresponding 
canonical momentum component, = pp^^ = —ihd/d(j) is the ^-component of 

the canonical orbital angular momentum, e = — |e| is the electron’s charge, and m its 
mass. 

We are interested in the evolution of non-stationary states of the system, described 
by the time-dependent Schrodinger equation 

( 2 ) 

where = (p, (j)) is the position of the electron in the plane perpendicular to the 
magnetic field. Note that if the electron is moving along the z axis with a velocity Vz^ 
the Schrodinger equation (2) describes the state of the electron after a propagation 
distance of 2 ; = Vzt [27]. 
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Figure 1: Cyclotron trajectory of the centre of mass of the wavefunction. This orbit 
occurs with respect to the axis y = yo, the position of which depends on the initial 
transverse momentum {p^^)o as well as the magnetic field. Also indicated is the 
direction of the canonical orbital angular momentum {Ll^^)z. 


An electron with momentum transverse to the magnetic field will exhibit cyclotron 
motion. This is conventionally described in a classical context. Here we will derive 
the cyclotron motion by assuming an electron wavefunction 


T(t = 0) = To = u{p) exp 




(3) 


where £ e Z. We have defined the x axis as the direction of the transverse kinetic 
momentum at t = 0. This state has a rotationally symmetric probability density 
|Top = |i4(p)p, and an expectation value of canonical angular momentum = 

£h. We shall see that the momentum = Pc^ results in a cyclotron orbit of the 

wavefunction. 

Note that in our model the canonical orbital angular momentum is not collinear 
with the instantaneous direction of propagation of the wavefunction. This is illus¬ 
trated in figure 1. The angular momentum is in the direction of the magnetic field, 
while the kinetic momentum has a component perpendicular to the magnetic field. 
This contrasts with vortex states, either in field-free space or in a magnetic field, 
which are energy eigenfunctions, as these have momentum and angular momentum 
which are collinear [29, 31]. Electrons in non-stationary states can have angular mo¬ 
mentum at an arbitrary angle to their direction of propagation, however [22, 32, 33]. 
Here, as a result of the cyclotron orbit, the time-averaged expectation value of kinetic 
momentum is collinear with the angular momentum 


3 Electron trajectories and angular momentum 

In the following we will show that the different forms of angular momentum give rise 
to rotations with respect to more than one axis. This can be seen by examining 
the “trajectories” associated with the electron’s probability distribution and current 
density. 

First, we will consider the expectation value of the electron’s position, which is 
equivalent to the centre of mass of its probability distribution. Differentiating twice 
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with respect to time, we obtain the equation of motion 


dHr±){t) 

dt‘^ 


-ul (rx -[x+ p^“/(eB)] x - [y - p^‘“/(eB)] y) (t) 
-u}^i{r±){t)-yoy), (4) 


where ujc = —eBjm is the cyclotron angular velocity and = Pc/{\e\B). Here we 
have used the fact that the quantity [x x ^ [y — which is 

the centre of the orbit of a classical particle which has position r± and momentum 
^kin constant expectation value yoy. The initial position and velocity of 

the centre of mass of the probability distribution are given by {r±){0) = 0 and 


respectively, and substituting these into (4) yields the trajectory 
{r±){t) = 2 / 0 [sinWei® + (1 - cosL 0 ct)y]. 


(5) 

(6) 


This trajectory, illustrated in figure 1, is a circular orbit with radius 


o- = |yo| = 


eB 


(7) 


and angular velocity uJc - the cyclotron orbit of a classical particle with the momentum 
Pc. The trajectory of the centre of mass of the probability distribution is therefore 
independent of the canonical angular momentum of the electron. 

The canonical angular momentum is instead associated with a circulation of cur¬ 
rent within the electron’s probability distribution. This can be seen by examining the 
probability current density j±(r_L, t) = Re(T*p^^^T)/m. To do so we have solved the 
time-dependent Schrodinger equation (2) numerically using the Chebyshev method 
[35, 36, 37], as described in Appendix A. We must first specify the radial distribution, 
u{p)^ of the initial wavefunction (3). Here we will set this to be the same as that of 
a Landau state - one of the energy eigenstates of the system: 


Hp) = ^n,\e\ 



(8) 


where n = 0,1,2,... is the radial quantum number, Ln is an associated Laguerre 
polynomial, pB = ^/4:h/\eB\ is the width of the Gaussian envelope and A^n,|^| = 
A/2n!/[7r(n-h \i\y.]/pB is a normalisation constant. This means that if Pc = 0, the 
electron would be in a Landau state. An arbitrary radial distribution could be de¬ 
composed in terms of the eigenfunctions . 

The time-evolution of the probability density |Tp and the current density j± are 
shown in figure 2. Here the transverse momentum pc has been chosen such that the 
radius of the cyclotron orbit is approximately equal to the width of the probability 
distribution. In (a) and (b) the electron has no net canonical orbital angular momen¬ 
tum, while in (c) and (d) it has a canonical orbital angular momentum i = 1. The 
evolution of these states is shown for different directions of the magnetic field, which 
result in different directions of the cyclotron orbit. Whereas the probability density 
follows a straightforward classical orbit, the current density is seen to depend in a 
non-trivial manner on both the wavefunction and the magnetic field. In particular. 
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in contrast to the classical cyclotron trajectory, and also to orbital angular momen¬ 
tum eigenstates in the absence of a magnetic field, the current distribution here is 
not rotationally symmetric. The rotational symmetry of the probability distribution, 
with respect to its centre of mass, is preserved, however. This reflects the fact that 
the magnetic field is rotationally symmetric, and the canonical angular momentum 
is conserved. This means that the canonical orbital angular momentum of the 
electron is associated with a rotation axis at the centre of mass of the probability 
distribution, and is independent of the cyclotron orbit. 

More surprising, perhaps, is that rotations in fact occur with respect to the cy¬ 
clotron axis and the centre of mass axis even when the electron does not possess 
any net canonical angular momentum. This is as a result of the diamagnetic angular 
momentum that any electron wavefunction possesses in the presence of a magnetic 
field [23, 24]. This angular momentum arises as a result of the circulating current 
the magnetic field induces within the wavefunction, and is associated with a rotation 
of the probability density at the Larmor angular velocity cjl = uJc/‘^ = —aBj(2vn) 
[27, 38]. The diamagnetic angular momentum is equal to = /'cjl, where 


r = m{p'^) 


(9) 


is the moment of inertia of the electron’s probability distribution, in the reference 
frame of its centre of mass [25]. Here p' = \r± — (r^)! is the radial coordinate in 
this reference frame. This angular momentum has the same direction as the external 
magnetic field, meaning that the associated magnetic moment, = eI/^^^/(2m), 
opposes the external field, in accordance with Lenz’s law. In contrast to cyclotron 
motion, the diamagnetic rotation depends on there being an extended probability 
distribution, and vanishes in the classical limit. For the Landau radial distribution 
the mean square radius is = (2n + \i\ + l)p^/2, and the diamagnetic 

angular momentum therefore takes the quantized values = sign(H)(2n+ \i\ + 

l)h [29]. The effect of the diamagnetic angular momentum becomes clear when we 
consider a superposition of opposite values of canonical orbital angular momentum, 
such as that shown in figure 3. As this superposition has no net canonical angular 
momentum, the rotation of the electron’s probability density with respect to its centre 
of mass is due entirely to the diamagnetic angular momentum. Previously we have 
interpreted this as a form of Faraday rotation for electrons [30] . 

In general, the orbital angular momentum with respect to the centre of mass axis 
will be given by a sum of canonical and diamagnetic contributions. The motion is 
thus described by two independent rotations: the cyclotron orbit, and the rotation 
around the instantaneous centre of mass axis due to the canonical and diamagnetic 
angular momenta. 

4 Parallel axis theorem 

The rotation of the electron’s wavefunction is reminiscent of a classical rigid body. 
We will explore this analogy further by considering the kinetic angular momentum 
of the electron, which is the total mechanical angular momentum it possesses while 
moving in the magnetic field. This angular momentum has the z component 



( 10 ) 


with an expectation value of 


{Lf^) = m + iujL 


( 11 ) 
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Figure 2: Probability density and current density j± for wavefunctions with 
canonical orbital angular momentum ^ = 0 (a, b) and ^ = 1 (c, d), for opposite 
directions of the magnetic field. In each case the transverse momentum is Pc = ^hjpB 
and the wavefunction has the radial distribution The red arcs indicate the 

trajectory of the centre of mass of the probability distribution, which is highlighted 
in green. 
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Figure 3: Probability density and current density for an equally weighted superposi¬ 
tion of wavefunctions which both have transverse momentum Pc = 2h/pB^ but have 
opposite values of canonical angular momentum ^ = dzl. The constituent states have 
the same radial distributions as in figure 2. 


for any state with a canonical orbital angular momentum Here I = m{p^) is the 
moment of inertia of the electron’s probability distribution for rotation with respect 
to the z axis. Just as with a rigid body, we can use the parallel axis theorem to 
express the moment of inertia / as a sum of two components: 

I = mPo + I’, ( 12 ) 

where po = |(t’±)| = \/2(l — coswct) a, with a defined by (7), is the radial coordinate 
of the centre of mass, and I' is the moment of inertia with respect to the centre of 
mass axis, given by (9). These two components correspond to the cyclotron orbit 
of the wavefunction and its diamagnetic angular momentum respectively. The total 
kinetic angular momentum of the electron, which we obtain from (11) and (12), can 
therefore be expressed as 


=m + Lf^ (13) 

where = muji^pQ = { 1 —cos L0ct)mL0c(T^ is the angular momentum associated with 

the cyclotron orbit. While the relation between the kinetic and canonical angular 
momenta in (10) is true also for a classical point particle [39], the decomposition 
into separate cyclotron and diamagnetic components which follows from (12) is only 
meaningful for an extended probability distribution. 

Unlike the canonical and diamagnetic components, the cyclotron angular momen¬ 
tum depends on our choice of reference axis. A natural choice is to consider the 
angular momentum with respect to the centre of the cyclotron motion. In this ref¬ 
erence frame, which we reach by making the transformation y y = y — yo, the 
cyclotron angular momentum has the constant value 

= mwco-^. (14) 

Irrespective of the reference frame, of course, the total kinetic angular momentum 
will be equal to the sum of the three components described. 

One may expect that the cyclotron angular momentum, which exists classically, 
would be the dominant contribution to the electron’s kinetic angular momentum. 
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However, this need not be the case, as can be seen by considering typical parameters 
for free electrons in electron microscopes. For example, if an electron beam which is 
initially propagating parallel to a magnetic field is transmitted through a diffraction 
grating with a period d = 100 nm, the first diffraction order will have a net transverse 
momentum of Pc = = hjd^ which corresponds to an energy of p\l(2m) = 

Ir ?= 0.15 meV. In a magnetic field of H = 1 T, the resulting cyclotron orbit, 
which has radius 41 nm, will have an angular momentum of 2.6h. This is of the same 
order of magnitude as the canonical angular momentum of the lowest order vortex 
states, and considerably smaller than that of vortex beams recently generated with 
a winding number of = 200 [21]. Indeed, in a given magnetic field, the cyclotron 
angular momentum can in principle have any size, ranging from zero to macroscopic 
values, depending on the net transverse momentum Pc. The diamagnetic angular 
momentum can also take a wide range of values, as the mean square radius of the 
probability distribution is varied [25], although this has a lower limit due to the 
uncertainty principle and a maximum due to the requirement of spatial coherence. 
This means that the different rotations we have described can indeed occur on the 
same length scale, justifying the choices of parameters in our figures. 


5 Superposition of cyclotron orbits 


So far we have considered rotations with respect to two different axes - the cyclotron 
axis, as well as the centre of mass axis. The position of the cyclotron axis was defined 
by the transverse momentum Pc which appears in a plane wave factor in the wave- 
function (3). Suppose, however, that we have a superposition of different transverse 
momenta. Arbitrary distributions of transverse momentum could be created using 
appropriately designed holograms [40, 41, 42, 43]. A simple example would be the 
following superposition of two momenta, Pc,i and Pc, 2 , which may also be associated 
with different canonical orbital angular momenta, £i and ^ 2 - 


To = 


V2 


exp 





+ 


V2 


exp 



Pc,2^ 


(15) 


This is similar to the wavefunction formed when a plane wave is transmitted through 
a forked diffraction grating [10]. For the discussion here the particular form of the 
radial distribution is not important, and the Landau function has been chosen with 
numerical efficiency in mind. The evolution of such a state, here with equal and 
opposite values of both the transverse momentum and canonical angular momentum, 
is shown in figure 4. It can be seen that there are now two cyclotron orbits, which have 
different rotation axes. These are associated with different directions of the initial 
momentum pc. As a result, the two components of the superposition move apart, 
before re-combining and interfering. 

Taken together, the cyclotron orbits in figure 4 describe a rotation, with respect 
to the 2 ; axis, at the Larmor angular velocity cjl- This is consistent with the predic¬ 
tions of classical electron optics regarding image formation in rotationally symmetric 
magnetic lenses [3]. Interestingly, though, in our case the axis of the Larmor rota¬ 
tion is not defined by a symmetry of the magnetic field - a uniform magnetic field is 
rotationally symmetric with respect to an infinite number of axes. Rather, here the 
Larmor rotation occurs with respect to the centre of mass of the electron’s probability 
distribution. This is the case both in figure 3, where the centre of mass follows a cy¬ 
clotron orbit, and in figure 4, where the centre of mass is stationary. Further, it must 
be remembered that we are considering here a single electron which is in a state of 
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Figure 4: Probability density and current density for an equally weighted superpo¬ 
sition of a wavefunction with Pc = 3h/pB and £ = —1 and a wavefunction with 
Pc = —3h/pB and £ = 1. The constituent states have the same radial distributions 
as in figures 2 and 3. The centre of mass of the superposition remains stationary 
in the x-y plane, while the individual components follow the cyclotron trajectories 
indicated. 


superposition. This means, for example, that if one of the two cyclotron components 
in figure 4 underwent an interaction which modified its phase, this could be detected 
through its effect on the subsequent interference pattern. 

6 Summary and outlook 

In summary, we have shown that in a magnetic field an electron can rotate around 
more than one axis simultaneously. The wavefunction of the electron follows a cy¬ 
clotron orbit, and superposed onto this is a rotation around the instantaneous centre 
of mass. The rotation with respect to the centre of mass axis arises as a result of the 
diamagnetic angular momentum, as well as any canonical orbital angular momentum 
the electron possesses. The kinetic angular momentum of the electron is therefore 
described by the parallel axis theorem. 

Our results show that canonical orbital angular momentum and cyclotron motion 
provide separate degrees of freedom for shaping electron current distributions. This 
could allow electron beams to be structured for use in specific applications. For 
example, the symmetry of the current distribution could be optimized to probe specific 
transitions in materials [40, 44]. It may also be possible to utilise cyclotron trajectories 
in novel forms of interferometry. Moreover, here we have only considered the case in 
which the canonical angular momentum and magnetic field are parallel, so that the 
rotation is confined to a plane. With canonical angular momentum and magnetic 
fields which are in different directions to one another, the angular momentum and 
current density could be shaped in three dimensions. 

Further, if the angular momenta were in different directions, it appears that they 
would become coupled. Canonical orbital angular momentum which is at an angle 
to a uniform magnetic field would be expected to precess around the direction of 
the field [22, 33]. The canonical angular momentum is also not conserved when the 
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rotational symmetry of the magnetic field is broken, such as in astigmatic magnetic 
lenses [45, 46, 47]. Not only this, but in non-uniform magnetic fields the spin and 
orbital degrees of freedom of an electron with non-relativistic velocity are no longer 
independent [12, 48]. The nature of the coupling between all of these angular momenta 
is an interesting avenue for future investigation. 
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A Numerical solution of Schrodinger equation 

The time-dependent Schrodinger equation can be solved numerically with high accu¬ 
racy and efficiency by expanding the time-evolution operator in a series of Ghebyshev 
polynomials. In this appendix, we describe how this method can be applied to the 
two-dimensional Schrodinger equation for an electron interacting with an external 
magnetic field. Eor generality, we shall consider here a Hamiltonian of the form 

«(».!/) = + (“) 

d 

+iSiy{x,y) — + So(x,y) 

[37]. In the case of the Hamiltonian used in the main text, we would have S 2 = 
—Ji?l{2m)^ Six = —heBy/{2m)^ Siy = heBx/{2m) and So = e^B‘^{x^ +^^)/(8m). 
The Mathematica code we have used to perform these calculations has been made 
available online at [50]. 

Since the Hamiltonian (16) is independent of time, we can write the solution of 
the Schrodinger equation as 


T(t + At) = exp ^(0- (17) 

In order to evaluate this numerically, first we must represent the wavefunction, and 
the coefficients S 2 etc., on a two-dimensional grid. If this grid covers an area LxX Ly, 
and contains Nx x Ny points, then the maximum spatial frequencies represented are 
kx,ma^ = 7 ^Nx/Lx and /c^^max = '^^y/Ly. The maximum and minimum values of 
energy represented on the grid are then 

-E'max ~ ^2(^a^^max T ^y,max) T ^^^x(^Six)kx^mayi (1^) 

+Max(5'i^)/c^,niax + Max(S'o) 

and 

^min = - [Max{Six)kx,ma^ + Max(5'i^)/c^,niax] + Min(5'o)- (19) 

Let us now introduce a new operator 
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where a = (^max ~ ^min)/2 and b = (^max + ^min)/2. This operator has eigenvalues 
represented on the grid which he in the range [—1,1]. We can then expand the time- 
evolution operator in a series of Chebyshev polynomials Tq{'H): 


T(t + At) = exp(—-6At) exp(—-a?iCAt)T(t) 

M 

^ exp(—^6At) (21) 

g=0 


Chebyshev polynomials are chosen as these minimise the error associated with trun¬ 
cating the expansion at a finite order M [35]. The expansion coefficients are given 
by 


aq{aAt) = 


{-i)^Jq{aAt/h), q = 0 
2{-iy Jq{a At/h), 


( 22 ) 


where Jq is a Bessel function. For q > eaAt/{2h), where e, in roman font, denotes 
the mathematical constant, the magnitudes of these coefficients decay exponentially 
with increasing q [49, 51]. This means that the error due to truncating the series at 
order M, which can be estimated by \aM{ciAt)\^ can be made arbitrarily small. If we 
set 

M = -^^ciAt (5, (23) 

6 can be adjusted so that this error is less than machine precision. The numerical 
error resulting from the Chebyshev expansion is then negligible. 

In order to evaluate the individual terms in the expansion (21), the action of the 
Chebyshev polynomial Tqij-L) on the initial wavefunction T(t) must be calculated. 
Using the recurrence relation for the Chebyshev polynomials, we obtain 


Tq{n)^{t) = 2nTq_,{n)^{t) - Tq_ 2 {n)^{t). (24) 


for g' > 0, with the initial conditions To(7i()T(t) = T(t) and Ti(?i()T(t) = ?iT(t). 
The action of the Hamiltonian on the wavefunction can be efficiently calculated by 
evaluating the spatial derivatives in Fourier space [49, 52]. That is, 

n^{t) « S2FT-^[{-kl - fc2)FT«'] + i5i^FT-i(ifc^FT^') 

+iSiyFT-\ikyFT'^) + So^, (25) 


where kx^ ky are the coordinates in Fourier space and FT denotes a discrete Fourier 
transform and FT“^ the corresponding inverse transform. 
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